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The tendencies to phase-separation and stripe formation of the t-J model on planes and four-leg 
ladders have been here reexamined including hole hopping terms t' , t" beyond nearest-neighbor sites. 
The motivation for this study is the growing evidence that such terms are needed for a quantitative 
description of the cuprates. Using a variety of computational techniques it is concluded that the 
stripe tendencies considerably weaken when experimentally realistic t' < 0, t" > for hole-doped 
cuprates are considered. However, a small t' > actually enhances the stripe formation. 



PACS numbers: 74.20.Mn, 74.25.Dw 

Growing experimental evidence suggests the exis- 
tence of static stripe order in a variety of transi- 
tion metal oxides, including hole-doped La2Ni04 and 
La1.6-xNdo.4SrxCuO4.OJ In other cuprates, such as 
Lai.85Sro.i5Cu04 and YBa2Cu30g.6) the magnetic scat- 
tering is consistent with the presence ofj-dynamic an- 
tiphase antiferromagnetic (AF) domains a The origin 
of these charge inhomogeneities is controversial. Some 
authors believe they are caused by Jahn- Teller (JT) 
distortions!!! Others favor a purely electronic explana- 
tion. For instance, hole domain walls were observed in 
Hubbard model Hartree-Fock calculations .□ In addition, 
computational calculations for the two-dimensional (2D) 
t-J model have found striped tendencies consistent with 
neutron scattering experiments.H'El The existence of diag- 
onal domain walls in four-leg t-J ladders has also been 
reported,0 adding to the expected strong similarities be- 
tween ladders and planes .El A third possibility is based on 
frustrated phase separation (PS) where the stripes arise 
from a combination of a slaprt-range attraction and long- 
range Coulomb repulsion. □ 

In parallel to these developments recent angle-resolved 
photoemission (ARPES) experiments addressed the one- 
particle spectral function of the undoped insulator 
Sr 2 Cu02Cl2.li3 The overall bandwidth and features along 
the (0, 0)-(tt, 7t) direction were found in agreement with 
theoretical t-J model predictions. However, the re- 
sults along (0, 7r)-(7r, 0) were puzzling since the ARPES 
quasiparticle-like peak has a clear energy maximum at 
(tt/2, 7t/2), while in the t-J model this line is almost flat. 
This difference is important and needs to be addressed. 

The main explanation proposed for the t- J- ARPES 
discrepancy is based on the relevance of corrections in 
the form of electronic hopping amplitudes beyond the 
nearest-neighbor (NN) contribution. The importance of 
these terms was recognized from the analysis of the elec- 
tronic structures of cupratesO In all these calculations it 



was concluded that for a proper description of cuprates 
a next-NN (NNN) hopping oLstrength t' along the pla- 
quette diagonal was necessary.E 2 ] For hole-doped cuprates 
t' has been systematically found to be of negative sign 
in contrast to the NN-hopping amplitude t with positive 
sign, and of about 20% to 40%-ii.s magnitude. Electron- 
doped cuprates need t' > 0x3 Note that in regions 
where AF correlations are important, t is renormalized 
to smaller values while t' is not severely affected, thus 
enhancing the relevance of such NNN corrections.. In 
addition, soon after the Sr 2 Cu0 2 Cl2 ARPES dataIL3 be- 
came available it was reported by Nazarenko et alx3 that 
including a t' < NNN-hopping the agreement theory- 
experiment was noticeably improved at all momeata. A 
similar result was also found by Lee and ShihJl3 Fur- 
ther work confirmed and improved this initial approach, 
showing that with the addition of an extra NNN i"-term 
connecting sites at distance,-2,a (a= lattice spacing) the 
results improved even moreilS Actually with these NNN 
corrections the hole-density dependence of the t-J spec- 
tra oLsmall clusters was found in agreement with ARPES 
dataO In Ref. 14 it was remarked that there is no sym- 
metry argument favoring the special case t' = 0, as it oc- 
curs in gauge theories where local symmetries and renor- 
malizability arguments fix the Hamiltonian. Since the 
existence of corrections to the t-J model are natural, the 
assumption t' = is mainly aesthetical. 

The goal of the present paper is to address the much 
discussed tendency of the t-J model to phase-separate 
and/or form stripes performing the calculations in the 
presence of a realistic nonzero NNN-hopping amplitude. 
The t-J Hamiltonian employed here is defined as 

H tJ = ~* E ( C 'W + H.c.) + J^(Si • Sj - mnj/4), 

<ij><7 (ij) 

where (ij) are NN sites. No doubly occupancy is allowed, 
t is defined as positive, and the rest of the notation is 
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Jit 

FIG. 1. (a) Stable and unstable regions of the extended t-J 
model on the in)- J/t plane. The results are obtained exactly 
on a 4 x 4 cluster with "averaged" boundary conditions (see 
text); (b) binding energy vs J/t of the extended t-J model on 
a 4 x 4 cluster, averaged boundary conditions, and using the 
NNN hoppings indicated; (c) site labeling used for the 20-site 
cluster. 

standard. The contribution of the NNN terms is 

H VV >=-t' ^ cl a C ma -t" ^ C t<r c r* + H.C., 
«km))<r («nr)»<7 

where ((km)) denote a pair of sites along the diagonals 
of the elementary plaquettes, and (((nr))) are pairs of 
sites located at distance 2a along the main axis. 

To find information about the influence of t' and t" 
on PS and stripe formation on the t-J model, here a va- 
riety of computational techniques have been employed. 
Let us start the analysis using the exact diagonalization 
(ED) method on^lanar systems. On 2D t-J clusters 
evidence in favorEa and againsto PS at small J/t has 
been recently presented. Our goal is not to add to this 
discussion, but rather to follow some of the approaches 
proposed in those papers and find the influence of NNN 
terms on the results. Let us start using ground state (GS) 
energies calculated by averaging with equal weight over 
a large number of twisted boundary conditions to reduce 
size efiects£3 Using this method the GS energies corre- 
sponding to the electronic density (n) = 1, the density 
under investigation (with N e electrons), and the density 
corresponding to N e + 2 electrons were obtained for in- 
creasing values of J/t. When the three energies lie on a 
single line, PS occurs between half-filling and the average 
between the densities N e /N and (iV e + 2)/NH 

Figure 1 (a) illustrates the influence of a NNN hopping 
on the PS tendencies of the t-J model using a 4 x 4 clus- 
ter. The largest density shown here corresponds to 13/16 
and was obtained with information from 12, 14, and 16 
electrons. For t' = the line separating the stable and 
unstable regions-converge to a very small J/t as the den- 
sity reaches 1.0. Eil However, even for an apparent "small" 
NNN hopping such as t' = —0.2t, the PS line now con- 



FIG. 2. Gijki for the seven configurations with the largest 
weight in the ground state (see text) vs J/t. (a) Results for 
t' — 0.0; (b) same as (a) but for t' = — 0.2t; (c) same as (a) 
but for t' = +0.2t; (d) same as (a) but for t' = +0.6t 

verges towards a larger J/t at half-filling. The effect is 
similar for t' = — OAt. Although these small-cluster re- 
sults should be considered only qualitative, the tenden- 
cies observed are clear and in agreement with a variety of 
calculations reported below. Then, it is apparent that a 
nonzero t' < amplitude reduces the tendencies towards 
PS in the t-J model near half-filling and small J/t. 

Figure 1(b) contains the (n) — 1 binding energy de- 
fined as A B = E(2)+E(0)-2E(1), where E(n) is the GS 
energy (averaged over boundary conditions) for n holes. 
The results show that a t' < term reduces the attraction 
between holes. This is correlated with a reduction of the 
probability of stripe formation and PS [Fig. 1(a)] caused 
by the short-range AF-induced attraction. The pairing 
region isj-pjushed up in couplings by a t' < since it fol- 
lows PS.lj Results qualitatively similar to those shown 
in Figs. 1(a) and 1(b) have been obtained using 18 and 
20-sites clusters and other values of t'r3 

To search for the dominant GS-hole configurations the 
matrix element G ijk i = {0\n h (i)n h (j)n h (k)n h (l)\0) was 
calculated. Here (i,j,k,l) denote four sites of the 20- 
site cluster following the labeling convention shown in 
Fig. 1(c). The product of the four hole- number oper- 
ators nh is a projector that provides the GS weight of 
the configurations with the holes at (i,j,k,l). The con- 
figurations with the largest Gijki are a — (4, 7, 14, 17) 
(unbounded holes), b = (3,7,13,17) (pairs of diago- 
nally bounded holes), c = (6,7,13,14) (NN-hole pairs), 
d = (6,7,13,18) (another NN-hole pair configuration), 
e — (2,8,13,19) [four-hole stripe along the (1,1) direc- 
tion], / = (1,2,3,6) (another type of four-hole domain 
hole), and g — (3,8,12,17) [four-hole stripe along the 
(1,0) direction]. Figure 2(a) contains the results in the 
absence of NNN hoppings. In this case three regimes can 
be identified: (1) at small J/t the holes are unbounded; 
(2) for J/t ~ 0.5 holes form pairs; and (3) at J/t above 
0.7 "diagonal" stripes are preferred, as observed first in 
Ref. 5. Although these results are qualitative, the ten- 



dencies towards (l,l)-stripe formation are clear. and also 
in agreement with four-leg ladder calculations.tj 

Figure 2(b) contains the results found for the same hole 
configurations but now using t' = —0.2t. Once again, in 
spite of the naively "small" value of t' its influence on GS 
properties is important. The (1,1) stripes are no longer 
competing with hole pairs and unbounded holes. Now 
the most relevant stripe configuration is the (1,0) stripe 
which dominates only for J/t ~ 1.9 or larger. Results 
similar to those shown in Fig. 2(b) have been obtained 
using a variety of t' < and t" > amplitudes. Thus, it 
is clear that the stability of the stripes is sensitive to the 
presence of NNN-hopping amplitudes. Since such hop- 
pings are expected to be realistic, the presence of stripes 
in electronic models for the cuprates with short-range 
interactions is called into question. 

It is interesting to note that using t' > 0, i.e., the 
"wrong" sign for hole-doped cuprates but relevant for 
electron-doped cuprates,E3 the tendencies to stripe for- 
mation are actually enhanced roughly in the small win- 
dow < t'/t < 0.2 at all values of J/t. Now the 
crossing point between configurations a and / appears 
in Fig. 2(c) at J/t <~ 0.4, while in Fig. 2(a) it occurred 
at J/t ~ 0.7. Then, a small and positive t'/t can be 
used as a test ground of electronic models with tendencies 
towards stripe formation. As t'/t grows further stripes 
become unstable again, and actually for t'/t around 0.5 
or larger the configurations with NN-hole pairs dominate 
for all the values of J/t explored here [Fig. 2(d)]. 

To understand the different influence of the sign of t' 
on the t-J model phase diagram, a discussion on the sub- 
tleties involving bare versus renormalized parameters is 
needed. It is known that at t' = 0, the effective NN- 
hopping amplitude is dramatically reduced at half-filling 
since intersublattice hole hopping distorts the AF back- 
ground. In this same regime effective nonzero t and 
t" amplitudes ace generated, as deduced from the one- 
hole dispersion £-3 The sign of this effective t' which gives 
mobility to the dressed hole is negative, a well-known 
fact which manifest itself in the minimum of the hole- 
quasiparticle band at (ir/2,ir/2). Adding a bare t' term 
of the same sign will enhance the hole mobility substan- 
tially, to the point where hole superstructures become 
unstable. This amplification of effects explains the re- 
sults of Figs. 1(a) and 2(a) and 2(b). However, the ad- 
dition of a bare t' in the Hamiltonian of positive sign 
can lead to a cancellation of effects, and a concomitant 
reduction of the hole mobility near half-filling. This will 
indirectly favor hole clustering since such structures arise 
from a competition of the potential energy gained by AF 
attraction and the kinetic energy—Spin correlations are 
enhanced with respect to t' = 0jl3 since poorly mobile 
holes cannot scramble them. When t'/t is increased fur- 
ther eventually holes should become mobile again and 
the stripe tendencies will diminish. Figure 2(d) shows 
Gijki now for t'/t = +0.6. The dominant configurations 
have NN-hole pairs instead of stripes at all values of J / 1 
shown. 
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FIG. 3. (a) e(x) (see text) vs x for several powers in the 
PL method, two values of J/t, and using a 6 x 6 cluster with 
t' — —0.2t; (b) same as (a) but using a 8 x 8 cluster. 



The GS energy of the t-J model supplemented by a 
NNN term with t' — — 0.2t has also itf en calculated us- 
ing the Power Lanczos (PL) method.t3 To search for in- 
dications of PS, the approach used in Refs. 19 and 26 is 
applied, namely the energy of the PS state is written as 
E = N s eii+Nhe(x), where N s is the number of sites, Nh 
the number of holes, eu the energy per site of the Heisen- 
berg model, x the hole density in a hole-rich phase, and 
e(x) — [eh(x) — en]/x. If at a fixed coupling J/t, e(x) 
is found to have a minimum at some density x m and the 
overall density is smaller than x m , then the system phase 
separates between a hole-free phase and a hole-rich one 
with density x m . In Fig. 3 e(x) is plotted vs x using 
6x6 and 8x8 clusters with periodic boundary condi- 
tions (open-shell configurations). The energies denoted 
by PL0-V correspond to results using an optimized trial 
wave function taken from the set of Gutzpdller and reso- 
nant valence bond (RVB) wave functions. E3Ej PL1-V de- 
note improved results now using the first Lanczos ste p ap- 
plied to the previously optimized wave function.EZI PL1- 
Pn correspond to further improvements resulting from 
the application of n powers of the Hamiltonian over the 
PL1-V wave function. For additional details the reader 
is referred to Refs. 19 and 25. 

At J/t = 1.0 and after the application of six pow- 
ers, the minimum of e(x) in Fig. 3(a) is found to be 
at x = 0.22. At J/t = 0.6, the minimum shifts with 
increasing powers towards the smallest density studied 
here namely x — 0.056. Similar results were obtained 
for the intermediate coupling J/t = 0.8 (not shown). In 
addition, using an 8 x 8 cluster analogous trends are ob- 
served as shown in Fig. 3(b), although with larger error 
bars due to the sign problem. The minimal values of e(x) 
are at the lowest doping density for J/t = 0.6 and 0.8. 
Therefore it is concluded that the critical J c /t for PS in 
the low hole density limit is at least « 0.8, which is larger 
than in the t' = case where J c /t rj 0.6.t3 The trends 
found here are qualitatively the same as observed using 
ED [Fig. 1(a)], namely a t' < moves the PS region 
towards larger J/t's. 
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FIG. 4. Rung-density (n r ) vs r on a 4 x 14 cluster, 
J/t = 0.5, 8 holes, studied with DMRG and open bound- 
ary conditions. Shown are only half the rungs, the rest is 
found by reflection, r = 1 (r = 7) is the end (middle) of the 
cluster. Results for several t''s are shown. 



Since the most robust computational evidence of stripe 
formation in one-band electronic models actually comes 
from the density-matrix numerical renormalization group 
(DMRG) studies of four-leg t-J ladders, let us complete 
our analysis by studying the same clusters as in Ref. 7, 
following a similar methodology, but now adding NNN 
terms. The DMRG results reported here were obtained 
using J/t = 0.5, m = 500 states, and a truncation error 
~ 1 x 10~ 4 . In Fig. 4 the rung density (n r ) (i.e., the sum 
of the four-site densities forming a rung) is shown using a 
4 x 14 cluster and 8 holes. Results for just half the lattice 
are provided for simplicity, since the rest are obtained by 
reflection. At t' — 0, previous resultsQ were reproduced 
as a test. They present a single broad peak indicative of 
GS clustering tendencies. Analyzing the hole-hole cor- 
relations the (1, 14 stripes were found to have a large 
weight in the GSE However, when the NNN amplitude 
is turned on, the effect is weaken. Consider, for example, 
t' = —0.2t: now a two-peak structure is observed which is 
more suggestive of hole pairing than of stripe formation 
(note there are four holes in average on the portion of the 
cluster shown in Fig. 4). The effect is further enhanced 
for t' — — Q.3t where the two peaks are sharper. Here the 
holes were found to be unbounded residing in pairs at the 
extremes of the same rung. Then, as grows a rapid 
transition from stripes to unbounded holes is observed. 
The melting of the stripes (assumed to be signaled by 
the melting of the single-peak in (n r )) roughly occurs at 
t' ~ — O.lt. More realistic values of the coupling are diffi- 
cult to study as accurately as at J/t = 0.5. However, it is 
expected that the short-range AF attraction will become 
weaker as J/t is reduced, and the tendencies to stripe 
formation will also be weaker in this regime. Note that 
once again a t > maintains the stripe structure (sin- 
gle peak shown in Fig. 4 for t' — +0.2t), as in planar 
systems. 

Summarizing, using a variety of computational tech- 



niques it has been shown that the tendencies to phase 
separation and stripe formation previously reported 
on planes and four-leg ladders are actually substan- 
tially weaken once realistic NNN-hopping amplitudes are 
added to the t-J model. The reason is that t' < terms 
enhance appreciably the mobility of holes, melting hole 
superstructures. Reciprocally, using a (small) t' > 
the stripes become more stable providing an interesting 
model for the analysis of their properties. 
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